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ABSTRACT 


To ensure sufficient capacity to handle unexpected demands for electric power, decision 
makers often over-estimate expeditionary power requirements. Therefore, we often use 
limited resources inefficiently by purchasing more generators and investing in more 
renewable energy sources than needed to run power systems on the battlefield. 

Improvement of the efficiency of expeditionary power units requires better manag¬ 
ing of load requirements on the power grids and, where possible, shifting those loads to 
a more economical time of day. We analyze the performance of a previously developed 
optimization model for scheduling time-shiftable electrical loads in an expeditionary power 
grids model in two experiments. One experiment uses model data similar to the original 
baseline data, in which expected demand and expected renewable production remain 
constant throughout the day. The second experiment introduces unscheduled demand and 
realistic fiuctuations in the power production and the demand distributions data that more 
closely refiect actual data. 

Our major findings show energy grid power production composition affects which 
uncertain factor(s) infiuence fuel consumption, and uncertainty in the energy grid system 
does not always increase fuel consumption by a large amount. We also discover that the 
generators running the most do not always have the best load factor on the grid, even when 
optimally scheduled. 
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Executive Summary 


Operational energy demand and requirements are major factors to consider in any military 
operation. The demand continues to grow as the U.S. military addresses conflicts across 
the globe and cooperates with other nations on humanitarian aid efforts. In fiscal year (FY) 
2013, the Department of Defense (DOD) used more than 9 million barrels of liquid fuel to 
support Operation Enduring Freedom in Afghanistan [1]. To ensure an uninterrupted flow 
of electrical energy for these military missions, most of the expeditionary and other DOD 
forward operating bases (FOBs) rely on fossil fuel-based generators to power equipment 
and systems for military operations. Over estimated power requirements account for this 
high-energy demand, which often results in an inefficient power grid system. 

Naval Postgraduate School Operations Research student John G. Sprague provides us with 
a tool to effectively manage load requirements in his thesis [2]. However, his model treats 
unmanaged load requirements and production as deterministic parameters. This thesis aims 
to investigate the impact of uncertain factors influencing fuel consumption and to obtain 
useful insights on the model’s performance to support future development of a robust version 
of it. 

We conduct a multifactor-designed experiment of key input parameters in Sprague’s model 
to more fully understand contributions to fuel savings by running two experiments. In 
both experiments, we use a nearly orthogonal Fatin hypercube (NOFH) design with 33 
design points. The baseline for the first experiment is very similar to Sprague’s baseline 
model but introduces systematic variation into six factors related to uncertainty in the 
system. In the second experiment, we seek additional insights into the model’s behavior 
in situations where the renewable energy is a larger part of the energy production and the 
unscheduled demand varies during the day. We change the composition of the power grid 
by increasing the renewable energy production and also modify our expected unscheduled 
demand and expected renewable production data to reflect patterns that are more realistic 
in an expeditionary power grid system. We conduct our analysis of these two experiments 
using JMP Pro Version 12.0.1 [3]. 




Our analysis pinpoints the heat transfer intereept as the dominant uneertain faetor eon- 
tributing to the eumulative total fuel eonsumption. The heat transfer intereept represents 
the shelter’s equilibrium heat transfer rate, whieh determines how fast heat is added to 
or removed from the shelter. Insights obtained from our study inelude the effeet of the 
eomposition of the power grid’s energy produetion on generator fuel eonsumption and the 
importanee of properly matehing the FOBs’ power requirements with the generators’ power 
ratings. Some notable findings inelude: 

• In terms of fuel eonsumption, the first experiment illustrates that the rolling horizon 
with perfeet future knowledge (RH-PFK) model variant always performs better than 
the rolling horizon with uneertainty (RH-U) by up to 8%. Also, the seeond experiment 
demonstrates that the RH-PFK performs better than the RH-U by up to 14%. 

• It is possible to meet power demand requirements with fewer resourees, sueh as the 
number of generators running simultaneously, and henee, less fuel eonsumption. 

• The seleetion of the type and eharaeteristies of shelters or struetures to ereet at the 
FOBs eould play a key role in fuel eonsumption. 

• The generator that runs most often does not always have the best load faetor. The 
mixture of generator ratings available for power produetion on the power grid affeets 
generator load faetors and subsequent fuel eonsumption. 

This study eomes from a desire to improve U.S. military energy use on the battlefield. 
The eost assoeiated with too mueh energy use goes beyond monetary resourees; lives are 
lost, and the future eapabilities budget is diminished. This thesis provides useful insights 
on how to reduee power grid energy eonsumption by using a powerful eombination of an 
optimization model and design of experiment (DOE) analysis. Ultimately, the U.S. military 
expeditionary power grid system should be resilient and self-sustaining for our assigned 
missions around the world. 
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CHAPTER 1: 

Introduction 


The demand for energy continues to grow as the U.S. military must address conflicts across 
the globe and assist with humanitarian assistance missions. Operational sustainability 
requires optimal energy use, especially in rugged and extreme environmental conditions. 
The increasing energy cost is a burden on the defense budget, taking away resources 
that could be invested in improving our military capabilities. Additionally, uncertainty 
associated with environmental conditions influences the design of an effective expeditionary 
power grid. The allocation of limited military resources to establish forward operating 
bases (FOBs) requires careful analysis of energy needs and resources. 

This thesis builds on the mathematical optimization model in [1] which exploits physical 
characteristics of heating ventilation and air conditioning (HVAC) systems and the energy 
sources that power them. Our study evaluates uncertain factors contributing to temperature 
differences in an expeditionary environment, using a design of experiment (DOE) approach 
to gain insight into the performance of the optimization model under a variety of circum¬ 
stances. Some of the findings may also support the development of a robust version of the 
optimization model. 

1.1 Department of Defense and Energy 

The world’s largest energy consumer, the Department of Defense (DOD), uses more energy 
in daily operations than any other private or public organization [2]. Comprising is not 
the right word here. The DOD consumes more energy than all other federal government 
agencies combined, which amounts to approximately 80% of the total U.S. government 
consumption [3]. The DOD’s energy use falls into two distinct categories: operational 
and installation or facility energy. In fiscal year (FY) 2014, operational energy accounted 
for 70% of the department’s energy usage; DOD used over 87 million barrels of fuel, at a 
cost of nearly $14 billion [4]. Figure 1.1 shows the DOD’s energy consumption over time. 
In 2008, the DOD accounted for 0.8% of the nation’s energy consumption, which is the 
equivalent of burning 395,000 barrels of oil per day. To put that figure into perspective, the 
DOD consumed about as much energy as needed to power the entire country of Greece [2]. 
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Figure 1.1. U.S. Government Energy Consumption by Agency (Defense vs. 
Non-Defense), Source [5]. 


To ensure the uninterrupted flow of electrical energy for military missions, most of the 
expeditionary and other DOD FOBs rely on fossil fuel-based generators to power equipment 
and systems for military operations. Although renewable energy currently flows into the 
power grid, fossil fuel-based energy far outstrips all other sources for DOD operational 
energy. Over a decade of war in Afghanistan and Iraq has drawn greater attention to the 
need for effective management and use of operational energy. To limit fuel usage, the 
DOD is identifying ways to reduce consumption. Figure 1.2 depicts the breakdown of U.S. 
operational energy consumption for FY 2014 by domain, service, mission, and combatant 
command. 

Each service spends a substantial amount of its budget on electrical power generation for 
operations across the globe. Air operations represent most of DOD’s operational energy 
usage. The U.S. Central Command (CENTCOM) is also a major consumer of energy, given 
the Iraq and Afghanistan wars. Air Force operations represent more than half of the entire 
DOD operational energy requirements [6]. The Navy’s operational energy requirements 
come from air operations as well as ships sailing around the world to maintain the freedom 
of commerce at sea. In years past, the Army and the Marine Corps accounted for the 
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Figure 1.2. Operational Energy Use, FY 2014, Source [6]. 


lowest consumption annually among the branches of service; however, this changed after 
the two recent wars when fuel requirements to support electricity generation and fuel costs 
skyrocketed. In FY 2013, the department used more than 9 million barrels of liquid fuel 
to support Operation Enduring Freedom in Afghanistan [7]. Table 1.1 shows that as the 
operational tempo increases, the demand for fuel also increases. 


Category 

Peacetime 

Wartime 


OPTEMPO 

OPTEMPO 

Combat 

Vehicles 

30 

162 

Combat Aircraft 

140 

307 

Tactical Vehicles 

44 

173 

Generators 

26 

357 

Non-Tactical 

51 

51 

Total 

291 

1040 


Table 1.1. Army Fuel Consumption in Peacetime and Wartime (Million 
Gallons per Year), Source [8]. 
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Table 1.2 displays some interesting observations of researehers in the Expeditionary En¬ 
ergy Assessment Environmental Control Unit Alternatives Study conducted for the Marine 
Corps Expeditionary Energy Office. The Marine Air-Ground Task Eorce (MAGTE) Power 
and Energy Model, (MPEM) uses data from Afghanistan and the operational cost of fuel 
at $6.39 per gallon to calculate savings in dollars. The result demonstrates that a 10% 
increase in efficiency from the current environmental control unit (ECU) suite requires only 
a nominal investment and the gain is technically feasible [9]. The subsequent increase 
in efficiency produces substantial savings that cannot be ignored. However, such savings 
naturally depend on the unit’s mission and the theater of operation. The 2009 field data 
from the USMC Afghanistan study concluded that generators averaged a load of 30%; addi¬ 
tionally, the HVAC represents 75% of the electrical demand and wastes 50% on inefficient 
structure [9]. Optimizing the HVAC control system and improving the efficiency of living 
structures is one of the most effective ways to reduce energy consumption. We analyze a 
model optimizing HVAC control systems in this research. 



Annual Fuel Savings in Tankers and Dollars 


10% ECUs 

20% ECUs 

30% ECUs 

Solar ECUs 

Tankers of Fuel 

79 

158 

217 

487 

Savings (Millions) 

$2.42 

$4.85 

$6.65 

$14.95 


Table 1.2. Annual Fuel Savings in Tankers and Dollars, Source [9]. 


1.2 Cost of Environmental Temperature Control 

In military operations, the cost of environmental control sometimes exceeds the cost asso¬ 
ciated with the equipment used in the mission. Overall costs include the energy used to 
power the equipment as well as to sustain it. The sustainment costs in combat and hostile 
environments sometimes outweigh everything else when there are human lives involved. 
In Eeslie Hayward’s interview with the first-ever Assistant Secretary of Defense for Opera¬ 
tional Energy Plans and Programs, Sharon Burke, she stressed that the DOD is ultimately 
a business, and that energy is a variable cost that must be kept down [10]. However, unl ik e 
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business, the DOD often engages in warfighting in which energy is both a liability and an 
opportunity [10]. We treat this energy cost as a liability in terms of fuel consumption, and 
review the impact of environmental temperature conditions on fuel consumption in fulfilling 
military missions. 

1.2.1 Monetary Costs 

Among numerous other reasons, the volatility of the world oil market makes it highly 
desirable to reduce the demand for operational energy. In 2010, the DOD consumed nearly 
5 billion gallons of petroleum in military operations, costing $13.2 billion, a 255% increase 
over 1997 oil prices [11]. The volatile price makes it difficult to accurately estimate and 
budget for fuel costs [11]. Overall, every $10 increase in the price of a barrel of oil increases 
the price of DOD operations by $1.3 billion [2]. To put this value in context, each $10 price 
increase is equivalent to the loss of almost the entire procurement budget of the U.S. Marine 
Corps (USMC) [2]. These estimates also do not account for the cost of transporting fuel 
to desired locations, which varies depending on the modes of delivery and the situation on 
the ground. For example, an Office of Naval Research (ONR) report found that the fully 
accounted costs (including transportation, storage, and logistics expenses) for a gallon of 
gasoline on the battlefield in Iraq and Afghanistan actually ranged from $15 per gallon to as 
much as $400 per gallon [2]. In-flight delivery of fuel for the Air Force costs approximately 
$42 per gallon [2]. 

Additionally, electricity generation constitutes a significant portion of operational energy 
usage on the battlefield. A late 2009 study conducted by the Marine Energy Assessment 
Team (MEAT) in Afghanistan revealed almost 89,000 gallons of daily fuel usage; out of this, 
power generation accounted for about 32% of the consumption [12]. HVAC systems alone 
used approximately 75% of the power generated [12]. To put this amount into perspective, 
of about 28,500 gallons of daily usage, HVAC systems consumed 21,000 gallons, costing 
about $22 million annually at $2.80 per gallon for the Marine Expeditionary Teams at the 
EOB used in this study. Also, the study uncovered the drastic growth in cost per gallon on 
the battlefield when approaching the tactical edge. This makes our $2.80 per gallon quite 
an optimistic estimate. The cost growth reflects the logistics involved in delivery. Eigure 
1.3 displays this cost growth. 
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Figure 1.3. Cost Growth Approaching The Tactical Edge, Source [12]. 

This cost burden deprives the DOD of funds that eould be invested in other programs to 
enhanee eapabilities to deal with new threats to national seeurity. 


1.2.2 Other Costs 


The eost growth, when approaehing the taetieal edge also eomes with a human eost. The 
most frequently employed taetie during warfare is deny or disrupt resupply to the opposing 
forees. The weapon of ehoiee eommonly used on eoalition forees’ fuel eonvoys in Iraq 
and Afghanistan has been improvised explosive deviees (lEDs). This ehoiee makes human 
life another variable, depending on the level of sophistieation of the opposing forees, the 
terrain, and the weather of the battle spaee. DOD offieials reported that in June 2008 alone, 
a eombination of these faetors resulted in the loss of about 44 trueks and 220,000 gallons 
of fuel [13]. Figure 1.4 illustrates the human eost of fuel resupply in Iraq. 

This figure shows that, between July 2003 and May 2009, lEDs aeeounted for about 43% of 
U.S. fatalities, a reminder of the fuel resupply burden in warfare operations. The manpower 
and equipment requirements to proteet the logisties route represent another eost faetor. 
Minimizing the fuel requirements will free up these resourees for other missions. 
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■ U.S. Iraq fatalities by lED ■ lEO % U.S. total U.S. Iraq FataliDes 



Figure 1.4. U.S. Fatalities in Iraq by lED, Source [13]. 


1.3 Energy Technology Initiatives and Opportunities 

The DOD recognizes this opportunity for savings and is investing substantial resources in 
energy initiatives. As part of the FY2013 budget, the DOD requested more than $1.4 billion 
for operational energy initiatives [3]. The DOD’s five year (FY2013-FY2017) Future 
Year Defense Plan (FYDP) includes a total of about $8.6 billion for operational energy 
initiatives [3]. The aim is to reduce the demand for energy required to conduct operations 
around the world. Table 1.3 shows some of these DOD energy initiatives. 

The rapid fielding of equipment in theater is one initiative designed to streamline the 
deployment of equipment, such as more efficient generators, to reduce in-theater fuel 
consumption. The implementation comes with the challenge of generator load and efficiency 
management. 

Cyclical loads from operating HVAC systems at FOBs present a unique obstacle to achieving 
maximum efficiency of the installed energy system. Environmental conditions, equipment 
operations, and mission requirements often determine the load requirements that should 
be planned for when setting up the FOBs. As the load requirements vary throughout the 
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Initiative 

Description 

Establish a baseline of DOD 
operational energy consumption 

Gather reliable DOD-wide data on operational energy consumption, to serve 
as a foundation for analyzing DOD operational energy use and developing 
operational energy metrics. 

Defense Operational Energy Board 

This board oversees the execution of the Operational Energy Implementation 
Plan, including setting operational energy metrics. 

Rapid fielding of equipment in 
theatre 

Efforts to streamline the deployment of equipment (such as more-efficient 
generators) that reduces in-theatre fuel consumption. 

Operational Energy Capability 
Improvement Fund 

$19.5 million fund to spur technology innovations that reduce energy load at 
contingency bases, measure energy consumption in forward areas, and 
transform waste into energy. 

Energy Key Performance Parameter 

Establish a methodology for an Energy Key Performance Parameter (KPP) to 
be used in requirements gap analysis, requirements development, and 
acquisition programs. 

Fully burdened Cost of Fuel 

Provide the services with rton-binding guidance on the methodology and 
application of Fully Burdened Cost of Energy (FBCE) estimates as part of the 
life-cycle cost analysis for new capabilities during the acquisitions process. 

Energy in Procurement Contracts 

Work with DOD's office of Defense Procurement and Acquisition Policy to 
develop template language on energy performance for DOD contracts. 


Table 1.3. Major ASD (OEPP) Operational Energy Initiatives, Source [3]. 


day, effective administration of the unmanaged requirements plays an important role in the 
efficiency realized from the installed capacity. Improving the efficiency of the installed 
capacity is critical to minimizing the generators’ fuel consumption. 


1.3.1 Load Sizing Peak and Average 

Sprague describes a four-step process that expeditionary energy systems use when analyzing 
the load requirement for a typical FOB [1]: 

1. Determine the type, quantity, and characteristics of connected loads in each structure 
or facility [1]. 

2. Calculate the theoretical maximum power requirement of a facility if all connected 
equipment is simultaneously operating [1]. 

3. Adjust this maximum value by applying a demand factor (<1) representing an as¬ 
sessment of the realistic portion of facility equipment that would simultaneously 
operate [1]. 

4. Add allowances for future growth [1]. 
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The current system requires the demand for air conditioning to be weighted by a demand 
factor of 1.0 [ 1 ], requiring the energy system to meet the full power requirement at all times. 
Unfortunately, environmental weather conditions can be very unpredictable, often leaving 
the generator with excess capacity. This excess capacity can lead to inefficient expeditionary 
energy systems. 

Sprague’s generator output simulation revealed that, when connected loads operate at a 75% 
duty cycle, the average power demand is 21.6 kilowatt (kW) and the peak demand is 28.8 
kW, resulting in high generator usage. On the other hand, at a 25% duty cycle, the result 
is the opposite with substantial wasted capacity and the generator frequently running below 
a 30% load, leading to more fuel burned inefficiently [1]. Field data from Afghanistan in 
Figure 1.5 displays the historical generator load from electrical demand. 



TQG »1,60kW Avg Load Olympia, 120kw Avg Load TQG »2.60kW Avg Load - Total Camp Avg Load 

Figure 1.5. U.S. Forces Afghanistan Camp Electrical Demand Over a 96- 
Hour Period, Source [1]. 
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1.3.2 Generator Fuel Efficiency 

Fossil fuel generators eonvert the ehemieal energy of the fossil fuel into the eleetrieal energy 
used to power expeditionary energy systems. The generator effieieney is the ratio between 
the useful eleetrieity outputs in a speeifie time unit to the ehemieal energy supplied from 
the fuel within that same time. Often, this effieieney is less than 50% with most effieieney 
lost to heat generated during the eonversion proeess. However, with an aeeurate estimation 
of the load faetor, larger Taetieal Quiet Generators (TQGs) tend to be more effieient than 
smaller units. Figure 1.6 illustrates the relationship between load faetor and fuel effieieney 
by TQGs sizes [1]. 



Figure 1.6. U.S. Military TQGs Fuel EfFiciency as a Function of Load Factor 
and Generator Size, Source [1], 

Selecting the smallest generator capable of providing the required power achieves maximum 
efficiency [1]. Sizing generators with load requirements based on the peak power demand, 
rather than the average power demand, creates inefficiency during low- and average-demand 
periods, which happens regularly. 

1.4 Research Objectives 

This research study investigates the impact of uncertain factors influencing fuel consumption 
in an expeditionary power grid. We conduct a multifactor-designed experiment of key input 
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parameters in Sprague’s optimization model to tease apart speeifie eontributions to fuel 
savings. We vary the uneertain faetors to understand speeifie effeets on the eumulative 
total fuel eonsumption, as well as to obtain useful insights into generator behaviors in the 
model. Our analysis of the impact of uncertainty on Sprague’s optimization model may 
also support the future development of a robust version of the optimization model. 

1.4.1 Scope 

This research employs the principle of DOE using the nearly orthogonal Latin hypercube 
(NOLH) design to vary uncertain input parameters in Sprague’s model. We extend the 
research beyond Sprague’s baseline power grid set-up by changing the energy production 
configuration, and introducing patterns for renewable energy production and unscheduled 
demand. This allows us to analyze how well Sprague’s model handles varying, uncertain 
environmental conditions, as well as assess the results of operating under these conditions. 

1.4.2 Thesis Contribution and Outline 

This thesis complements industry and government progress in the field of efficient tacti¬ 
cal hybrid power systems and reduced energy demand. The analysis of optimized load 
scheduling provides more insight into practical considerations for using this optimized load 
scheduling approach. 

Chapter 2 presents a literature review of the academic, commercial, and government studies 
that inform our present research. In Chapter 3, we describe the model developed by Sprague, 
and our DOE methodology that investigates uncertain factors influencing generator fuel 
consumption. Chapter 4 presents our findings, and Chapter 5 contains the conclusion and 
future work. 
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CHAPTER 2: 
Literature Review 


Both the government and industry often diseuss the optimization of eleetrie power gener¬ 
ation from various energy sourees. The U.S. military, in partieular, has invested heavily 
in initiatives and researeh to improve effieieney on generators, ECUs, and struetures used 
at FOBs for eondueting military operations. In his thesis, Sprague highlights the intersee- 
tion of three fields of researeh and applieation: military energy efheieneies, hybrid smart 
mierogrids, and demand side management [1]. In eontrast, this thesis draws from the pre¬ 
vious work and evaluates uneertain faetors eontributing to temperature differenees in an 
expeditionary environment, using a DOE approaeh to gain insight into the performanee of 
the optimization model under a variety of eireumstanees. Some of the findings may also 
support the future development of a robust version of the model analyzed. 

We review a brief deseription of the expeditionary energy power system, energy effieieney 
trade-offs and power generation, distribution systems, and the prineiples of the design of 
experiments. 

2.1 Expeditionary Energy Power System 

A taetieal expeditionary energy power system will eonsist of some eombination of a power 
generation subsystem, a distribution subsystem, and the eleetrieal loads requirements of 
the eustomers. Also, energy storage is sometimes ineorporated into the power system 
to enhanee its resilienee and effieieney. Figure 2.1 depiets a general representation of a 
typieal expeditionary power system. Sprague provides a more detailed deseription of the 
expeditionary system and various energy system arehiteetures [1]. 

2.2 Energy Efficiency Trade-offs and Power Generation 

Reeent developments inelude the following: more effieient power generators; better insu¬ 
lated struetures; intelligent, software driven eontrollers for distribution; and more effieient 
demand side equipment drawing power from the grid. Energy effieieney trade-offs abound 
in eleetrie load management. Eleetrie load management is a speeifie way of eontrolling 
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1: Base tents 

2: Feeders and distribution centers 
3: Generators 


Figure 2.1. Expeditionary Energy Power System Set-Up, Adapted from [14]. 

the peak load in the network to produee a eonstant or average demand [15]. This process 
involves both the supply-side and demand-side management of power requirements. 

Efforts to reduce generator fuel usage and cost have motivated some academic research in 
recent years. A Naval Postgraduate School (NPS) student, R. Kelly [16] investigated the 
benefit of transitioning from traditional generator employment to an alternative architecture 
using an Energy Management System (EMS). The EMS provides an intelligent interface 
by matching real-time load demand to the smallest capable power source to reduce fuel 
consumption between various power sources, loads, and energy storage elements to form 
a microgrid [16]. T. Mansfield et al. [17] at the United States Military Academy created a 
prototype Hybrid Tactical Power System (HTPS) by drawing energy from various sources to 
address the low-peak load variability problem via load aggregation for a more efficient use 
of generators, therefore reducing fuel usage and generator maintenance cost. Bouaicha [18], 
a student at NPS, examined a Hybrid Electric Grid (HEG) model that uses weather forecasts 
to establish an optimal day-ahead schedule for the grid in a robust and cost efficient manner. 
Another NPS student, Ulmer [19], extended this research further with a capital planning 
optimization model by using historical solar data and simulated forecasts for wind data to 
formulate a mixed integer linear program that maximizes the islanding time of a microgrid, 
subject to budget and physical constraints. Islanding time is the amount of time demand 
can be met without a connection to the commercial power grid. 
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The commercial utility industry, the supply side, studied the demand pattern from con¬ 
sumers, the demand side, for over a decade to establish a pricing system that would influ¬ 
ence power demand requirements. Unmanaged peak load power demand requirements will 
eventually require utility company to increase generation capacity to meet the increasing 
demand. Increased capacity becomes excess during the minimum and average load periods, 
therefore creating an inefficient and expensive energy system. To reduce the peak load 
period and bring it closer to the average load period, the pricing mechanism incorporates 
an incentive pricing system. It requires the shiftable cyclic load to be delayed until the 
beginning of the minimum or average load demand period for better pricing. The design 
allows either the customers or the utility company to control load shifting during the peak 
period [15]. Figure 2.2 displays what the utility company aims to achieve with this approach. 
The end result is the demand side with a lower bill and the supply side with a reduced peak 
load, increased load factor, and improved efficiency [15]. 


Shifting loads from peak period between noon and 1600 around. 



0:00 4 00 8:00 1200 16:00 20 00 24:00 

Time of Day (Noun) 


Figure 2.2. Peak Load Shifting, Source [15] . 

On another note, efficiency trade-offs in the military have been mostly pursued with the 
introduction of more efficient generators and the addition of hybrid power plants [1]. From 
TQGs to Advanced Medium Mobile Power Sources (AMMPS) and Ground Renewable 
Expeditionary Energy Network Systems (GREENS), improving efficiency continues to 
dominate the DOD’s effort at reducing the energy demand at EOBs. The recent effort of the 
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U.S. Army “The Tactical and Deployed Power” includes the 1 kW multi-fuel generator, an 
alternative energy hybrid generator that provides intelligent, tactical, and convenient squad 
level power [20]. 

2.2.1 Renewable Energy Power Generation 

Renewable energy sources provide the expeditionary power generation system the resilience 
required to sustain longer operations without relying too much on fossil fuel energy. Al¬ 
though still being continuously researched, solar-powered unmanned aerial vehicles (UAVs) 
will provide the DOD a great option to reduce fuel consumption in expeditionary units. 
The top companies currently using this technology and their popular solar UAVs include 
the following: 

• The Airbus - Zephyr 

• Boeing Phantom Works - SolarEagle 

• Google Titan Aerospace - Solara 

• Facebook Ascenta - High-speed internet via drones 

• AeroVironment/NASA - Gossamer Penguin, Solar Challenger, NASA Pathfinder 
(Plus), NASA Centurion, NASA Helios 

• Lockheed Martin - Hale-D 

The Renewable Energy for Distributed Under-Supplied Command Environment (REDUCE) 
[20] is a program developed by the U.S. Army Communications Electronics Research, 
Development and Engineering Center. It is a state-of-the-art system enabling the harvesting 
and use of solar and wind power, energy storage, energy distribution and monitoring, and 
mobility; it can produce up to 5 kW of power in the field. As shown in Figure 2.3, it can 
be moved around easily with its light tactical flat deck trailer mount, making it easy to be 
placed in a position of maximum utility and performance. 

Solar-powered ECU, produced by SunDanzer, an El Paso-based commercial firm in business 
since 1999, employs solar panels to power a cooling only ECU when there is sufficient 
sunlight, without using batteries for storing excess power for future usage [9]. Figure 2.4 
shows a typical set-up for this system. The direct current air cooler (DCAC) model developed 
by National Aeronautics and Space Administration (NASA) and licensed to SunDanzer uses 
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Figure 2.3. REDUCE, Source [20]. 

a patented refrigeration teehnology in its design by optimizing the refrigeration eyele to 
aehieve low energy eonsumption [21]. 



Figure 2.4. Solar Powered ECU, Adapted from [9]. 

2.2.2 Other Power Generation 

The traditional approaeh to power generation uses fossil fuel-powered generators. The 
Projeet Manager for Mobile Eleetrie Power (PM MEP) is the DOD exeeutive agent that 
oversees the battlefield eleetrie power generation and integration. Eigure 2.5 shows the 
organizational ehart. PM MEP missions are as follows [22]: 
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• Establish, maintain, and standardize a DOD-wide family of Mobile Electric Power 
Generating Sources (MEPGS), including the development of the best mobile gener¬ 
ators. 

• Execute program, including consultation on how to improve and expand standardiza¬ 
tion MEPGS across the DOD. 

• Approve and disapprove all requests for non-standard MEPGS. 


Matnx Support 
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Figure 2.5. Team MEP Organization Chart, Source [22]. 


The PM MEP assisted in the DOD transition from the TQGs to AMMPS as one of the 
measures to improve effieieney in power generation. The Offiee of Naval Researeh (ONR) 
is spearheading eollaboration with aeademia, industry, and the PM MEP on the use of 
modular, eompaet solid oxide fuel eell (SOEC) generators for taetieal eleetrieal power, 
whieh will provide highly efiieient silent power generation that is towable for mobile taetieal 
units. 


2.2.3 Energy Storage and Power 

With an energy storage deviee built into the system, the power grid’s resilienee ean be 
improved. Reeent advanees in teehnology have improved the storage eapaeity and length of 
time for battery energy storage. The DOD eollaborates with industry, government ageneies. 
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and academia to improve energy storage for operational use. The Hybrid Energy Storage 
Module (HESM), shown in Eigure 2.6, is a collaborative effort between the DOD and the 
Department of Energy to develop a sealable storage platform for military and eommereial 
use [23]. More information on energy storage ean be obtained from Sprague [1]. 





Prime Power 
(Generator) 


Isolation 

Switch 


Loads 


Figure 2.6. Energy Storage Module, Source [23]. 


2.3 Distribution Systems 

The DOD eonducts operations mostly away from existing installations where infrastruc¬ 
tures already exists. The operation site’s energy requirements must be generated on-site, 
through an efficient power generation system, and distributed through an intelligent system. 
Eigure 2.7 illustrates sueh a distribution system. Previous design efforts foeused on how to 
supply and distribute power whenever required at the right quantity. Today’s distribution 
efforts address properly sizing and distributing energy to obtain an efficient energy system. 
It involves the design of an intelligent mierogrid system that is eapable of addressing the 
following: 

• Development of real time and non-real time models, simulation, and optimization to 
support designs that are adaptable to multiple eomponents, deviees, and systems on 
the grid [23]. 
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• Robust and flexible bi-directional power control at power densities suitable for mul¬ 
tiple platforms [23]. 

• Integrated hardware, open architecture design and intelligent automatic control, and 
a unified management system. 

• Power converters and inverters enabling interoperability. 



1. Gefi*r«tor Set 

2 . Feeder Center 

3. Distribution Center 

4. Field Tent 


Figure 2.7. Typical Distribution System, Source [24]. 


The opportunity for an intelligent power management system can be pursued through 
networking, advanced sensors, and control software and hardware. The future of the current 
Power Distribution Illumination System Electrical (PDISE) is intelligent power distribution, 
in which every component in the system is driven by model software that can sense and 
predict requirements on a continuous basis, and then distribute loads in the most efficient 
and optimal manner. Eigure 2.8 displays the ultimate goal. This thesis contributes toward 
this field of research. Our analysis of the impact of uncertainty on the optimal scheduling 
of an electrical load model will shed more light on the areas that require more attention. 
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Figure 2.8. Intelligent Power Distribution System, Source [14]. 


2.4 Design of Experiments 

Every day the world, ineluding systems and proeesses around us, eontinues to beeome more 
eomplex. Making any informed decision on designing future systems and processes requires 
understanding today’s complexities and the factors contributing to this complex behavior. 
Design of experiments (DOE) is a powerful tool that enables researchers to gain insight into 
cause-and-effect relationships that characterize our complex world. Central to the principle 
of DOE are models or systems with input parameters or factors, undergoing a process or 
transformation to produce one or more outcomes. Each outcome is called a response whose 
behavior depends, in part, on the input factors called the predictors. Scholars continue to 
investigate and research options for generating designs that allow analysts to gain insight 
efficiently and effectively when running experiments. 

DOE dates back to the work of Eisher [25], a mathematician whose pioneering effort 
mostly dealt with applications in agriculture. Some of his discoveries include the analysis 
of variance (ANOVA), principles of blocking and randomization, and Eatin square and 
split plot experiments [26]. Many designs intended for physical experiments focus, by 
necessity, on varying only a few factors and severely limiting the number of design points. 
More recently, large-scale designs suitable for large numbers of factors have proven to be 
very beneficial for simulation experiments; examples and references are provided in [27] 
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and [28]. For example, Sanehez et al. [27] highlight the importanee of using large-seale 
designed experiments in defense and homeland seeurity applieations. Their simulation 
experiment involving a hybrid optimization model for unmanned aerial vehiele surveillanee 
demonstrates this signifieanee, saving an estimated $20 billion for the U.S. Army [27]. 
Spaee-filling designs like the Latin hypereube, the orthogonal Latin hypereube, and the 
NOLH are useful. 

Applieation areas of DOE are expansive and eontinue to grow. As just one example, over 
160 theses and numerous papers involving simulation experiments appear on the web pages 
of the NFS’s SEED (Simulation Experiments & Efficient Designs) Center for Data Earming 
[29]. Those most closely related to this work evaluate optimization models. Eor example, 
Gardner [30] conducts a series of experiments on a two-stage stochastic optimization model 
for planning naval logistics functions in humanitarian assistance operations, and Ei [31] uses 
a similar approach to find a model for planning the military response following typhoons in 
Taiwan. Gardner and Ei both find that the results of changing the inputs to the optimization 
model in a designed experiment yields useful insights about its performance and robustness. 
Morse [32] investigates the robustness of an optimization model to predict future energy 
needs and costs for Marine Corps expeditionary operations, including the use of solar 
energy; the robustness of his results in the presence of correlated cost coefficients is further 
explored in Sanchez et al. [33]. 
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CHAPTER 3: 
Methodology 


Sprague’s “Optimal Scheduling of Time-Shiftable Electric Loads in Expeditionary Power 
Grids” model minimizes fuel consumption by prescribing an optimal schedule for generator 
operation, ECU utilization, and energy storage management in a deterministic manner. The 
model considers the physical constraints of generator and ECU operation, as well as thermal 
constraints that must be satisfied within the internal structures in expeditionary EOBs [1]. 

We first provide a brief description of Sprague’s model. We then lay out two designed 
experiments that assess the performance and robustness of this model under varying condi¬ 
tions. 


3.1 Sprague’s Optimization Model 

A detailed description of Sprague’s model can be found in Chapter 4 of his thesis, but we 
summarize its features to better explain our contribution to the research. 

Sprague’s model makes a set of assumptions for simplification; the thermal model is the 
principle constraint that must be observed to maintain the structure’s internal temperature 
within acceptable limits; and the formulation of the optimization model, which incorporates 
the assumptions and thermal model. 

3.1.1 Objective 

Sprague’s optimization model is a discrete-time mixed integer linear formulation that min¬ 
imizes generator fuel consumption. It uses prescriptive mechanisms for periods of ECU 
operation, subject to specific internal structure temperature requirements, and the physical 
limitations of equipment. 

The objective is to minimize total fuel consumption, subject to a number of constraints. The 
complete list of constraints with their general descriptions can be obtained from Sprague 
[1]. Some of the constraints include those for keeping the temperature within a specific 
range, those that ensure power requirements properly match power production, those that 
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limit generators’ start and stop times, and so on. Sprague studies several variants of the 
optimization algorithm, but we focus on two that seemed to perform well and reflect real-life 
application: rolling horizon with perfect future knowledge (RH-PFK) and rolling horizon 
with uncertainty (RH-U). Both models are run in a rolling horizon fashion, meaning that 
they iterate over time steps by optimizing over a planning horizon and then executing that 
plan over a shorter execution horizon. Following execution, the models then advance their 
planning horizon forward by the length of the execution horizon, repeating this process 
until all time steps have been considered [1]. The RH-PFK receives perfect knowledge 
during this process, meaning that the parameter values it encounters during the planning 
process are the same as it encounters during the execution process. The RH-U, on the 
other hand, receives a set of constant forecasted values for planning purposes and the actual 
values during execution. For demonstration purposes, Sprague generates actual demands 
from a specified uniform random distribution with constant expected value [1]. For our 
experiments, both models receive the same sets of actual parameter values. We describe 
our process for generating these values in Sections 3.3 and 3.4. 

3.1.2 Model Assumptions 

The following assumptions ensure model simplification and computational tractability and 
help to bridge the gap from shortfalls in data availability and computation power: 

• Loads are balanced in each period and commensurate with the total load of all 
operating generators split among them according to their nameplate rated capacity [ 1]. 

• The charge and discharge rate limits of energy storage devices are constant for the 
range of the battery level, and the level varies linearly with the chosen charge and 
discharge rates [1]. 

• All power production is aggregated and available to any load within each period [1]. 

• Software tools not independently validated with field or test data provide estimates of 
the thermal behavior of shelter and ECU performance [1]. 

• Generator fuel consumption is modeled as a piecewise linear curve of rated power [ 1]. 

3.1.3 Thermal Model 

At any point in time, the internal structure temperature is a function of the external environ¬ 
mental temperature or conditions, the physical structure characteristics, and the ECU work 
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output [1]. Although the model is similar regardless of whether heating or cooling is being 
considered, Sprague only considers cooling for the sake of simplicity [1], and we will do 
the same. 

3.1.4 Inclusion of Variability 

Sprague accounts for uncertainty variability with the introduction of three parameters: 
unmanaged demand, renewable production, and the heat transfer intercept. Both variants 
are run using the same sets of actual values, randomly generated from Uniform distributions 
centered around the expected values, with floors of zero to prevent negative results [1]. 
Sprague’s model parameters use the following uncertain data as inputs. 

• Unmanaged or unscheduled demand consists of demand from the customers that arises 
from the unplanned use of equipment such as ironing clothes, using the microwave to 
prepare meals, charging electronic devices, and so on. 

• Renewable production comes from sources that include solar, wind, hydroelectricity, 
and so on. Power production from these sources varies throughout the course of the 
day. 

• The heat transfer intercept represents the shelter equilibrium heat transfer rate, which 
determines how fast heat affects the shelter. The environment adds heat to the shelter, 
as do the equipment and personnel within the shelter. The intercept parameters 
depend on the characteristics and type of shelter. For instance, wooden and brick 
shelters will have different mean heat transfer intercepts. The intercept may vary 
during the course of the day as a result of people entering and leaving the shelter. 


3.2 Baseline Scenario 

Sprague uses data from Afghanistan to derive his model parameters [1]. We use many of 
the same values; our baseline scenario also uses a time step of two minutes, and assesses the 
performance over 300 time steps (a ten-hour day). This baseline sets up our first designed 
experiment in this research. 
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3.2.1 Model Baseline Composition 

The baseline grid eonfiguration models a U.S. Marine Corps (USMC) patrol base in 
Afghanistan with 45 Marines, sheltered in eight struetures, using 10 ECUs, three eon- 
ventional generators, and three hybridized generators with battery storage [1]. Figure 3.1 
provides a general overview of sueh a power grid. 
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Figure 3.1. Southwest Afghanistan USMC Patrol Base showing Equipment 
Inventory, Source [1], 


Sprague models this baseline grid with nine shelters, four generators, mean renewable 
energy produetion of 1.1 unit, and one battery storage with eharge and diseharge rate of 
0.2 times battery eapaeity, respeetively. The minimum allowable run and rest times for 
the generators are five minutes eaeh, and all generators and renewable power are available 
to serve any of the shelters’ eooling requirements. The aetual minimum run and rest 
times for the generators are rounded up to six minutes beeause they must be multiples of 
the time step. More details on the inventory eomposition, equipment, storage, renewables, 
unmanaged demand, and uneertainty ean be obtained from Chapter 5 of Sprague’s thesis [ 1]. 
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Tables 3.1 through 3.3 show the parameters for baseline shelter and generator eonfigura- 
tions, as well as fuel consumption. We use the same mix of shelters and generators as in [1], 
although the identifying numbers for the shelters and the generators are slightly different. 



iase Shelter Confignrations 

Parameter 

Shelters 1-4, 7 

Shelters 5, 6 

Shelter 8 

Shelter 9 

Usage 

Billeting 

Billeting 

Billeting 

COC 

T>pe 

Base-X 305 

Base-X 305 

Base-X 307 

Base-X 307 

Min. internal temp. [°F] 

70 

70 

70 

70 

Max. internal temp. [°F] 

80 

80 

80 

80 

ECU capacit>' [BTU hr] 

36,000 

20,000 

60,000 

96,000 

ECU power [kWT 

4.5 

1.6 

8.5 

13.0 

Thermal slope [BTL7(hr.°C)] 

-931 

-931 

-1,226 

-1,254 

Thermal intercept [BTU hr] 

37,790 

37,790 

48,794 

60,454 

ECU run / rest time [minutes] 

HI 

I'l 

2'2 

m 


Table 3.1. Model Baseline Shelter Configuration, Adapted from [1]. 


Baseline Generator Configurations 


Parameter 

Generator 1 

Generator 2 

Generator 3 

Generator 4 

iNpe 

-AMMPS 

-AMMPS 

-AMMPS 

-AMMPS 

Frequency [Hertz] 

60 

60 

60 

60 

Rating [Kw] 

10 

30 

60 

10 

Gen. run rest time [minutes] 

5/5 

5 ■' 5 

5 ■ 5 

5 ,'' 5 


Table 3.2. Model Baseline Generator Configuration, Adapted from [1]. 


FUEL CONSUMPTION USING AVL 


Power 

Percent 

Generator Model 

5-kW 

10-kW 

15-kW 

30-kW 

60-kW 

gal/hr 

gal/hr 

gal/hr 

gal/hr 

gal/hr 

110% 

0.55 

0.98 

1.39 

3.11 

5.33 

100% 

0.51 

0.88 

1.24 

2.79 

4.92 

75% 

0.42 

0.70 

0.95 

2.00 

3.96 

50% 

0.34 

0.53 

0.73 

1.39 

2.74 

25% 

0.27 

0.38 

0.49 

0.92 

1.66 

10% 

0.23 

0.29 

0.38 

0.65 

1.08 

0% 

0.20 

0.24 

0.31 

0.59 

0.74 


Table 3.3. AMMPS Generator Sets Fuel Consumption, Source [1]. 
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3.2.2 Measures of Effectiveness 

Our primary measure of effectiveness is cumulative fuel consumption. We examine this 
for both the RH-PFK and RH-U variants of the model, so another measure is the relative 
efficiency of RH-U to RH-PFK. We also look at the use of the generators, both in terms 
of total time, and the time operating at low efficiency. Additionally, we examine the 
temperature variation in each of the shelters in terms of the cooling performance of the 
ECUs. Both of our experiments assess these measures of effectiveness. 

The baseline results in Sprague’s thesis are illustrated using a single set of uncertain 
parameter data for his optimization model. We run multiple replications for a variety of 
different scenarios, chosen in structured ways, to provide new types of insight into the 
model’s performance and robustness. 

3.3 Experiment I: Baseline Scenario 

This experiment analyzes Sprague’s model and directly compares its results to his. It 
uses the same set of input parameters, including the number of shelters and the shelters’ 
parameters, the same number of generators, and so on. Uncertainty implementation is also 
similar in that it includes the renewable power production, the unscheduled demand and the 
mean heat transfer intercept; the only difference is in the heat transfer intercept deviation, 
which is lower. 

3.3.1 Factors 

Table 3.4 displays the input factors that bring uncertainty into the model. For each distri¬ 
bution, we use the three parameters to generate independent sample data for time step 0 
through time step 300. The “mean” represents the expected value. The “deviation” repre¬ 
sents the half-width of a uniform distribution, expressed in terms of the proportion relative 
to a baseline value. In experiment I, the baseline value is the mean for all three sources of 
uncertainty. 

Specifically, (3.1) shows the distribution of the unscheduled demand data at any given time 
step: 
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Unscheduled'^'^™""^ ~ (3.1) 

Unscheduled'”®"" * (l + Uniform [l - Unscheduled®^®", Unscheduled^®"]). 

We only consider Unscheduled*" < 1, so all generated values are greater than zero. 

The heat transfer intercept is also uniformly distributed, as shown in (3.2). 

Intercept"'""*" ~ Intercept'"®"" * (Uniform [l - Intercept*", 1 + Intercept*"]) . (3.2) 

Finally, the renewable production is uniformly distributed, as shown in (3.3). 

Renewable^""*®"'"" ~ (3.3) 

Renewable""®"" * (Uniform [l - Renewable*", 1 + Renewable*"]) 


3.3.2 Design 

We construct an NOLH design [34], using an Excel spreadsheet template avail¬ 
able at the Simulation Experiments and Efficient Designs (SEED) center website: 
https://harvest.nps.edu. Table 3.4 shows the uncertain factors and the ranges of values 
for each of the six factors. This design has 33 design points (i.e., combinations of factor 
settings), and can be used for up to eleven factors. The maximum pairwise correlation is 
0.242. We ran five replications of the design. The results appear in Chapter 4. 


Uncertain Factors 

Factor 

TMJe 

Ranges 

Unscheduled mean 

Continuous 

[6.7, 9] 

Unscheduled de\iation 

Continuous 

[0.1, 0.2] 

Renewable mean 

Continuous 

[0.9, 1.3] 

Renewable desiation 

Continuous 

[0.1, 0.2] 

Intercept mean 

Continuous 

[0.9, 1.1] 

Intercept de\iation 

Continuous 

[0.2, 0.5] 


Table 3.4. Factors for Baseline Experiment. 
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3.4 Experiment II: Time-Varying Production and Demand 

We add three additional faetors to ensure our unseheduled demand and renewable pro- 
duetion data refleet a more realistie pattern in an expeditionary power grid system. The 
faetors for the heat transfer intereept remain unehanged. However, instead of holding the 
mean unseheduled demand eonstant during the entire period, we use three faetors — the 
unseheduled base similar to unseheduled mean, the halfwidth used as the middle of the 
base, and offset used to eontrol the peak height — to introduee a pattern that shows de¬ 
mand rising and falling during the eourse of the day. Similarly, we use three faetors — a 
renewable magnitude multiplier, the halfwidth, and peak offset — to generate an expeeted 
mean renewable produetion that varies over the day. We use these faetors to introduee a bell 
shaped pattern of a single eyele of a eosine funetion that aeeounts for the periods of zero 
to little demand and produetion, and periods of peak or high demand and produetion. We 
superimpose uniform noise, where the noise deviation is a faetor (similar to experiment I). 
Figures 3.2 and 3.3 display examples of the trends and values we obtain. The trend in Figure 
3.2 represents a more realistie depietion of the power grid system in whieh power demand 
rises gradually in the morning, peaks during the afternoon when the temperature rises and 
ECUs must provide more eooling, and drops eventually at night. Renewable produetion 
follows a very similar pattern of high solar energy produetion during daytime, braeketed by 
zero produetion at night. 



Figure 3.2. Example of Non-constant Unscheduled Demand Pattern, Exper¬ 
iment II . 


To visually understand how this works, see the examples in Figures 3.3 and 3.2. We set the 
expeeted renewable energy to peak at 120, whieh represents the renewable offset, and set 
the halfwidth to 100, whieh means that the renewable drops down to zero at time 20 and 
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2.5 

2.0 

1.5 

1.0 

0.5 

0.0 


time 220, as shown in Figure 3.3. The aetual renewable produetion is set by superimposing 
a noise uniform random variation around the aetual renewable to reflect the inconsistencies 
in renewable production. The peak of the expected renewable production in the left plot is 
equal to 2.0, representing a case where the renewable magnitude multiplier is 1.0. A larger 
value of the multiplier will result in a larger expected renewable production. 

In Figure 3.2, we set the expected unscheduled demand at a halfwidth of 180, which 
ensures the unscheduled demand does not drop down to zero during the experiment. The 
actual unscheduled demand follows a similar pattern with a noise uniform random variation 
superimposed to reflect the unpredictability of unscheduled power demand. 

3.4.1 Factors 

The factor settings for this experiment are similar to those of experiment I, but now include 
the three new factors. Table 3.5 illustrates the factor settings. One other potential factor 
that we do not explore is the peak offset for the unscheduled demand curve; we keep it fixed 
at midday (time step 150), although this could be changed in future studies. 

The ranges for the unscheduled halfwidth ensure that the mean unscheduled demand never 
falls to zero during our ten-hour day. In contrast, for the renewable production we explore 
combinations of halfwidths and offsets that can fall to zero at either or both ends of the 
day, or remain positive throughout the day. Consequently, we choose to vary the renewable 
magnitude multiplier, halfwidth, and deviation of the mean as separate factors, and let these 
determine the total renewable production. We describe our process for using these factors 
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Figure 3.3. Example of Non-constant Renewable Energy Pattern, Experi¬ 
ment II. 
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Uncertain Factors 

Factor 

Type 

Ruees 

Remarks 

Unscheduled base 

Continuous 

[6-7, 9] 

Similar to die unscheduled mean 

Unscheduled halfwidth 

Integer 

[180, 300] 

Half\\'idth used to control the load peak 

Unscheduled de\iation noise 

Continuous 

[0.1, 0.2] 

Similar to the unscheduled de\Tation 

Renew'able magnitude 

Continuous 

[9, 13] 

Multiplier for the renew’able 

Renewable mean halfvvidth 

Integer 

[90, 210] 

Halftvidth used to control production drop-off times 

Renewable de\iation noise 

Continuous 

[0.1. 0.2] 

Similar to renewable de\iation 

Renewable peak offset 

Integer 

[120, 180] 

Also used to control production drop-off times 

Intercept mean 

Continuous 

[0.9, 1.1] 

Same as experiment I 

Intercept de\iation 

Continuous 

[0.01, 0.025] 

Same as experiment I 


Table 3.5. Factors for Secondary Experiment. 


as well as demonstrate pietorially in Seetion 3.4. Figures 3.4 and 3.5 show additional 
illustrations of how the faetors were used. 
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Expected unscheduled demand vs. time 



Setting the halfwidth at 
180 ensures 
unscheduled demand 
would not fall down to 
zero during the 300 
time steps considered. 
It reduced to zero at 
around *30 and 330. 
Increasing the 
halfwidth to 300 
spreads the curve out 
with the unscheduled 
falling to zero at 
around 150 and 450. 


Figure 3.4. Non-constant Unscheduled Demand Illustration, Experiment II. 


3.4.2 Design 

The NOLH design setup for the seeond experiment resembles the first experiment, exeept 
for the addition of the three new faetors. The maximum pairwise eorrelation is 0.0234. We 
ran five replieations of the experiments, and the analysis of the results appears in Chapter 4. 
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Setting the 
renewable offset 
to peak at 120 
and the 
renewable 
halfwidth set at 
100 means that 
the renewable 
drops down to 
zero at time steps 
20 and 220 to 
reflect times of 
day with no 
renewable 
production. 


Figure 3.5. Non-constant Renewable Energy Illustration, Experiment II. 


3.5 Summary 

We conduct two experiments using two 33-design point NOLHs to determine which of 
the uncertain factors most affect overall fuel consumption. These factors include the 
unscheduled demand requirements, the renewable energy production, and the heat transfer 
intercept. 
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CHAPTER 4: 
Analysis of Results 


We now present the results of our experiments, including the summary outcomes of 
Sprague’s model as well as the results from our DOE. Our baseline scenario with experi¬ 
ment I uses the same configuration as Sprague’s with the exception of a smaller intercept 
deviation factor value, so we compare the results obtained directly with Sprague’s in this 
DOE analysis. 

Also, we analyze the experimental data to gain additional insight into the behavior of 
the model, the uncertain input parameters, their relationships and interactions, and their 
influence on the cumulative total fuel consumption. 

4.1 Model Results 

As previously discussed, Sprague’s model minimizes fuel consumption by optimally 
scheduling which of the four generators and nine ECUs to operate, while maintaining 
the shelter’s internal temperature within the specified limit. The model results include 
the cumulative total fuel consumption, the generators’ operation schedule, and the ECUs 
operation schedule. 

4.1.1 Model Fuel Consumption Results 

Sprague uses the legacy system (EGCY) model variant as the baseline or the standard to 
measure model improvements. The EGCY allows the generators to run continuously without 
managing any resources, ensuring enough power to support all the connected loads [1]. It 
implements a simple threshold policy for the ECUs, turning on each ECU when its shelter 
hits an upper temperature limit and turning it off when the shelter’s temperature reaches 
a lower limit. Table 4.1 shows the cumulative fuel consumption of Sprague’s model for 
three model variants: the EGCY used as the baseline, the RH-PEK and the RH-U, over a 
10-hour optimization horizon. We also introduce a relative efficiency measure to compare 
the performance of the two model variants, RH-PEK and RH-U, that are the focus of this 
study. We compute this parameter as the ratio of cumulative total fuel for the RH-U to that 
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of the RH-PFK. Any time the relative effieieney is greater than one, the RH-PFK performs 
better. As expeeted, the results show that the LGCY model consumes the most fuel [1]. 
The RH-PFK and the RH-U record substantial improvements over the LGCY variant [1]. 


Baseline Cumulative Fuel Consumption - 10 hours 

Variant 

Best Solution 
[gallons] 

Best Possible 
[gallons] 

Reduction vs. LGCY 
[0/.] 

Relative EfBciencj' vs. 
RH-PFK Best Possible 

LGCY 

32.9 

32.9 

NA 

1.44 

RH-PFK 

22.9 

22.9 

30.4 

NA 

RH-U 

23.7 

23.5 

28 

1.02 


Table 4.1. Baseline Configuration Cumulative Fuel Consumption, Adapted 
from [1]. 


Design point 17 from our DOE experiment I closely resembles the single run in Sprague’s 
baseline. Not only is the baseline scenario the same, but so are the values of the unscheduled 
mean, the renewable mean, and the intercept mean. Specific values for the actual demand 
differ from those of Sprague, but generating the actual unscheduled demand and the actual 
renewable production requires using identical distributions. The only difference is that 
we use a narrower distribution for generating the actual intercept values: Sprague uses an 
intercept deviation of 0.05, while we use a value of 0.02. The values of a comparison of the 
cumulative fuel consumption results obtained from our experiment’s five replications are 
similar to those of Sprague. The values of RH-PFK range from 22.00 to 22.42, the values 
of RH-U range from 22.13 to 23.51, and the relative efficiency ranges from 0.55 to 2.50. 
Overall, the average cumulative fuel consumption of approximately 22.9 gallons for RH-U 
and 22.2 gallons for RH-PFK is only about 3.5% more fuel consumption for the RH-U as 
reported by Sprague. 


4.2 DOE Results 

The model results from General Algebraic Modeling System (GAMS) provide a substantial 
amount of information spread across multiple files. We use IMP Pro Version 12.0.1 to filter, 
consolidate, organize, and process our data into a format that can be easily analyzed [35]. 
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4.2.1 Initial Exploration of Experimental Data 

We first examine a subset of our data set at time step 300 to obtain cumulative fuel 
consumption at the end of the time steps for each model run and each design point. 

Figures 4.1 and 4.2 demonstrate that this subset of the relative efficiency distribution has 
outliers. Further investigation reveals that the outliers resulted from infeasible solutions 
obtained by the model run for some design points, although no design point yields infeasible 
solutions for more than two of its five replications. One potential explanation is that 
occasionally, the optimization solver “times out” on a particular iteration based on the 
tolerance time limit set for the model. This can lead to an infeasible solution for subsequent 
iterations of the rolling horizon algorithm. The time limit set for the model in each iteration 
of the rolling horizon algorithm is 3,600 seconds. The overall completion times for all 
iterations vary widely in our experiments, from 1,953 seconds to 10,757 seconds. 

If Sprague’s model is to be used in actual operations, then this infeasibility requires fur¬ 
ther investigation. A larger number of replications would provide more insight into the 
proportion of time the optimization model is unable to find a feasible solution in a timely 
manner. A major advantage of using DOE to study systems or models is its ability to vary 
factors simultaneously and gain insights not obtainable from varying such factors one at a 
time. This means that if certain factors or factor combinations contribute to the chances of 
obtaining feasible solutions, the DOE may help identify them. 

We exclude outliers from further consideration in our analysis, and the relative efficiency 
distributions without outliers appear in Eigures 4.3 and 4.4. As expected, RH-PEK consis¬ 
tently outperforms RH-U, by up to 8% in experiment I and up to 14% in experiment II, as 
shown in Eigures 4.3 and 4.4. On average, RH-PEK outperforms RH-U by 2.8% for ex¬ 
periment I. Therefore, the results are consistent with Sprague’s conclusion that uncertainty 
increases the fuel consumption only slightly. Our experiment validates that this perfor¬ 
mance holds whenever a solution is obtainable, and does so across multiple replications for 
all design points. Eor experiment II, RH-PEK outperforms RH-U, on average, by 5.1%. 

4.2.2 Generator Behavior Insights 

We are interested in the uncertain factors’ influence on the cumulative total fuel consumption 
as well as the behavior of other model responses. The generators’ operation and load factor 
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^ Rel Eff (Uto PFK) Distributions 



Figure 4.1. Experiment I Time 
With Outlier. 


Quantiles 

A 

Summary Statistics 

100.0% 

maximum 

7.487804S78 

Mean 

1.0673571 

99.5% 


7.487804878 

Std Dev 

0.5031477 

97.5% 


1 .0703948975 

Std Err Mean 

0.03917 

90.0% 


1.0557801835 

Upper 95% Mean 

1.1446996 

75.0% 

quantile 

1 .0367062407 

Lower 95% Mean 

0.9900146 

50.0% 

median 

1 .0270692348 

N 

165 

25.0% 

quantile 

1.0148458492 



10.0% 


1.0098141213 



2.5% 


1 .0059443863 



0.5% 


1 



0.0% 

minimum 

1 




Step 300 Relative EfFiciency Distribution 


Rel Eff (U to PFK) Distributions 



Quantiles A 

Summary Statistics 

100.0% 

maximum 14.196721311 

Mean 

1.1685921 

99.5% 

14.196721311 

Std Dev 

1.073992 

97.5% 

1.9990842491 

Std Err Mean 

0.0836102 

90.0% 

1.0876500639 

Upper 95% Mean 

1.3336832 

75.0% 

quartile 1.0658448579 

Lower 95% Mean 

1.0035009 

50.0% 

median 1.0484615385 

N 

165 

25.0% 

quartile 1.0352863726 



10.0% 

1.0247188145 



2.5% 

1.0129694961 



0.5% 

1.0096836669 



0.0% 

minimum 1.0096836669 




Figure 4.2. Experiment II Time Step 300 Relative EfFiciency Distribution 
With Outlier. 


Rel Eff (U to PFK) Distributions 



Quantiles A 

Summary Statistics 

1 00.0% 

maximum 1.07798871 

Mean 

1.028 2081 

99.5% 

1.07798871 

Std Dev 

0.016494 

97.5% 

1.06883104 

Std Err Mean 

0.001283 

90.0% 

1.05507383 

Upper 9 5% Mean 

1.0307513 

75.0% 

quartile 1.03649146 

Lower 9 5% Mean 

1.0256648 

50.0% 

median 1 .0270144 

N 

164 

25.0% 

quartile 1.01473109 



10.0% 

1.00979357 



2.5% 

1 .0059385 



0.5% 

1 



0.0% 

minimum 1 




Figure 4.3. Experiment I Time Step 300 Relative EfFiciency Distribution 
Without Outlier. 


affect overall fuel consumption, and a closer look at these may help to obtain useful insights 
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Figure 4.4. Experiment II Time 


Quantiles a Summary Statistics 


100.0% 

maximum 

1.1 324786325 

Mean 

1.0512081 

99.5% 


1.1324786325 

Std Dev 

0.0236721 

97.5% 


1.11 50037667 

Std Err Mean 

0.0018656 

90.0% 


1.0834770563 

Upper 95% Mean 

1.0548925 

75.0% 

quartile 

1 .0639902426 

Lower 95% Mean 

1 .0475237 

50.0% 

median 

1.0477941176 

N 

161 

25.0% 

quartile 

1 .0348962524 



10.0% 


1.0241274553 



2.5% 


1.0129634564 



0.5% 


1.0096836669 



0.0% 

minimum 

1.0096836669 




Step 300 Relative EfFiciency Distribution 


Without Outlier. 


into the generators’ behavior in the model. The GAMS output already ineludes eaeh 
generator’s status and load faetor during operation, for every time step. Each generator has 
two levels representing its status: zero when idle and one when running. Each generator’s 
load factor is also represented by a continuous variable between zero and one. Erom these 
variables, we can compute the number of generators operating at every time step. We can 
also obtain insights into the proportion of time each generator operates, the proportion of 
time that specific numbers of generators are turned on during the optimization, and the 
proportion of time generators operate below or above a load factor performance level. 


4.2.3 Experiment I 

Using the results of experiment I, we begin by analyzing the uncertain factors’ influence on 
the cumulative total fuel consumption using the partition tree method. Partition trees tease 
apart relationships between response variables and the predictor variables by recursively 
splitting the data set where relationships occur, thereby creating a tree-like structure. The 
heat transfer intercept mean multiplier of the thermal model is the only influencing factor 
after three splits for both the RH-PEK and the RH-U model variants, yielding RSquare 
values of 0.94 and 0.93, respectively. Eigures 4.5 and 4.6 display the results. RSquare is a 
statistical measure that helps to determine the proportion of variance in our response factor, 
or the cumulative total fuel, which can be explained or predicted by the predictor factors, or 
the uncertain factors. It is not surprising that the partition tree results are nearly identical for 
the two model variants, because the correlation between the two performance measures is 
nearly 100% (0.99). The influence of other uncertain factors is minimal in this experiment 
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since at least 93% of the variation eould be explained by the intereept mean alone. We 
eonclude that the heat transfer faetor (the intereept mean multiplier), has the most influence 
on the fuel eonsumption of the factors we tested and subjeet to the degree of variability. A 
higher heat transfer intereept means that more fuel is used. For instanee, if there are different 
types of shelters available, selecting those with lower heat transfer intercepts to reduee fuel 
eonsumption. Of eourse, operational eonsiderations may also constrain the choiees. 



Figure 4.5. Experiment I Factors’ Influence on Cumulative Total Fuel Parti¬ 
tion Tree, RH-PFK. 


Furthermore, we apply a linear regression approaeh to fit a model for the eumulative total 
fuel, by allowing stepwise regression to select among the main effeets, quadratic effects, 
and two-way interactions. This approach helps to investigate the type of relationship that 
exists between our response faetor and the predietor faetors. Figures 4.7 and 4.8 show the 
results. The outeomes parallel those of the partition tree approach for both model variants, 
but also reveal that the relationship between eaeh of the important faetors and fuel used 
is linear. The intereept mean dominates other faetors, followed by the unseheduled mean, 
and the renewable mean, in that order. Onee again, the intereept mean eould explain at 
least 92% of the variation in our response faetor. The side-by-side plots in Figures 4.7 and 
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Figure 4.6. Experiment I Factors’ Influence on Cumulative Total Fuel Parti¬ 
tion Tree, RH-U. 


4.8 highlight the difference in the RSquare value when we use one dominant factor as a 
predictor versus three significant factors. We conclude that the heat transfer intercept mean 
contributes the most to cumulative fuel consumption in the Sprague model’s approximation 
of the expeditionary power grid, ffowever, the majority of the power production capacity 
for the grid comes from the generator in Sprague’s model. Consequently, the renewable 
production would likely be more influential if it were represented by a larger proportion of 
the production capacity. 

Practically, decision makers need to invest in better shelter structures with reduced heat 
absorption properties, as well as develop control measures that prevent heat transfer to the 
shelters to reduce fuel consumption. Measures such as controlling cooling boundaries that 
ensure individuals transiting shelters close the doors plays an important role in reducing 
fuel consumption. None of the other uncertainty factors has any major influence on the fuel 
consumption because the majority of the power generation is from the generator, and (as we 
shall discuss shortly) the unscheduled demand is not putting any major stress on the power 
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grid. Thus, the optimization model ean handle that amount of uneertainty explored in the 
DOE whenever it eould find a feasible solution. 


Effect Summary: Three Factors' Predictors 


Effect Summary: One Factor Predictor 


Source 

intcpt_mean_mult 

unsched.mean 

renew.mean 


LogWorth 


198,481 

102.515 

14.614 



PValue 

0.00000 

0,00000 

0.00000 


Remove Add Edit Undo I I FOR 


Source LogWorth 

intcpt_mean_mult 99.2781 
Remove Add Edit Undo I I FDR 


Actual by Predicted Plot 


Actual by Predicted Plot 



RSq=1.00 RMSE=0.1327 


Summary of Fit 

RSquare 

0.996749 

RSquare Adj 

0.996688 

Root Mean Square Error 

0.132674 

Mean of Response 

22,16634 

Observations (or Sum Wgts) 

164 



RSq=0.94 RMSEsO.5736 


Summary of Fit 

RSquare 

0.938477 

RSquare Adj 

0.938098 

Root Mean Square Error 

0.573577 

Mean of Response 

22.16634 

Observations (or Sum Wgts) 

164 


PVatue 

0.00000 


Figure 4.7. Experiment I Factors’ Influence on Cumulative Total Fuel Step¬ 
wise Regression, RH-PFK. 


Additionally, we examine the generators’ behavior in the model to obtain some insights 
on use and performanee. Table 4.2 shows the proportion of time eaeh generator operates. 
Both model variants use the 10 kW generators, generators 1 and 4, about the same amount 
of time (about 34% for RH-PFK and about 33% for RH-U). The range of use for both 
generators is also about the same. It appears that the model uses smaller eapaeity generators 
to handle loads first, then turns on the larger generators if the smaller ones eannot handle 
the load. This is expeeted sinee larger generators are less efiieient when operating at low 
loads. We eonelude that the model performs well in managing fuel efheieney by matehing 
load requirements with smaller generators first and then using larger ones when the load 
requirement exeeeds their eapaeity. Generator 2 is the least used for both model variants at 
about 17%. 

Also, the proportion of the number of generators turned on, as shown in Table 4.3, suggests 
that the model rarely requires more than two generators at any point. A single generator 
operates approximately 69% of the time for RH-PFK and approximately 70% of the time for 
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Effect Summary: Three Factors' Predictors 


Effect Summary: One Factor Predictor 
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Figure 4.8. Experiment I Factors’ Influence on Cumulative Total Fuel Step¬ 
wise Linear Regression, RH-U. 


Proportion of Time Generator Operates 


RH-PFK 

RH-U 

Generator 

Proportion in % 

Range 

Proportion in % 

Range 

1 (lOkWO 

34.76 

23.59 - 48.17 

33.07 

17.28 -46.84 

2 (30k\\0 

17.44 

7.64 - 32.89 

17.58 

8.64-31.23 

3 (60k\\0 

39.35 

24.92 - 51.50 

47.17 

29.90 - 61.13 

4 (lOkWO 

34.30 

22.59 - 50.50 

32.92 

22.59 -48.17 


Table 4.2. Experiment I Proportion of Time Each Generator Operates, RH- 
PFK vs. RH-U. 


RH-U. The results display that the power grid’s need for four generators at any point in time 
is almost zero. Further studies need to investigate whether the power grid can efficiently 
operate with fewer generators and still be fuel efficient. 

Lastly, we analyze the load factor performance of the generators as shown in Figures 4.9 and 
4.10. The quantile summary statistic reveals that 75% of the operating loads of the small 
generators (generators 1 and 4) are at least an 85% load factor for the RH-PFK. This result 
indicates that, when used, these two generators operate at a high load factor which reduces 
fuel consumption. The situation is similar for the RH-U when 50% of generators 1 and 4’s 
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Proportion Number of Generators Turned On 


RH-PFK 

RH-U 

No of Generator 

Proportion in % 

Range 

Proportion in % 

Range 

0 

3.24 

0.33-8.97 

0.92 

0.00 - 4.65 

1 

68.93 

55.81-79.07 

69.63 

56.48 - 83.06 

2 

26.60 

15.61-40.53 

27.29 

16.28 - 39.20 

3 

1.20 

0.00-4.32 

2.13 

0.00 - 7.97 

4 

0.03 

0.00-0.99 

0.04 

0.00 - 0.99 


Table 4.3. Experiment I Proportion of The Number of Generators Turned-On 
Every Time-Step, RH-PFK vs. RH-U 


loads are at least an 81% load faetor. This difFerenee exhibits that generator 3’s load faetor 
is not as good as generators 1 and 4’s load factors, particularly for RH-U. We conclude that 
a further experiment could vary the number and capabilities of the available generators in 
order to see whether different configurations result in better load factors, while still meeting 
the demand. Realistically, the installation of the right combination of generator ratings 
improves power grid load factor efficiency, which directly translates into fuel savings. 
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Figure 4.9. Experiment I Distribution of Generators’ Load during Operation, 
RH-PFK. 
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Figure 4.10. Experiment I Distribution of Generators’ Load during Operation, 
RH-U. 
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4.2.4 Experiment II 

This experiment introduces more realistic unscheduled demands and renewable production 
data into our analysis. Also, the experiment incorporates a larger amount of renewable en¬ 
ergy production into the grid. We employ the same evaluation approach as for experiment I. 

Once again, in a few cases the model does not report a feasible solution for the RH-PFK 
(specifically, one replication for design points 2 and 13, and two replications for design 
point 6). We omit the infeasible replications in the analyses that follow. We analyze the 
uncertain factors’ influence on the cumulative total fuel consumption using a partition tree. 
The expected intercept and the renewable halfwidth with the RSquare value of 0.80 are the 
only influencing factors after three splits for the RH-PFK. However, for the RH-U model 
variant, the influencing factors are the renewable halfwidth, the expected intercept, and the 
renewable noise deviation with an RSquare value of 0.66. Figures 4.11 and 4.12 display 
the results. The difference in outcomes from the experiment I results signifies that the grid 
configuration, demand patterns, and renewable production pattern have a definite impact on 
fuel consumption. Additionally, we see that fuel consumption savings result if we increase 
renewable energy production on the power grid. 

In addition, we use stepwise linear regression to fit a model for the cumulative total fuel at 
time step 300, allowing for main effects, quadratic effects, and two-way interactions. This 
approach helps to uncover the type of relationship that exists between our predictor factors 
and this response. Figures 4.13 and 4.14 show that simple relationships perform well for both 
of the optimization model variants. For both model variants, the influencing factors include 
the expected intercept, the renewable halfwidth, renewable magnitude multiplier, and the 
unscheduled base. RH-PFK recorded a RSquare value of 0.99 and RH-U also documented a 
RSquare value of 0.99, an indication that almost 100% of the variation in the cumulative total 
fuel could be attributed to these influencing factors. This experiment II outcome highlights 
the importance of the proportion of power production from various energy sources in grid 
power generation capacity. This outcome also shows that the uncertainty of the renewable 
production has greater influence than other sources of uncertainty for the factor ranges we 
investigate. The signs of the coefficients make sense. For instance, more fuel is used when 
the expected intercept increases while less fuel is consumed as we increase the renewable 
halfwidth. This implies that if we have a higher peak renewable production, we end up not 
using the generator as much which translates into reduced fuel consumption. 
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Figure 4.11. Experiment II Factors’ Influence on Cumulative Total Fuel 
Partition Tree, RH-PFK. 


Furthermore, we look into the generators’ behavior to obtain useful insights on performanee 
and utilization. Table 4.4 exhibits the proportion of time eaeh generator operates. Generators 
1 and 4, both with power ratings of 10 kW ran the most in the RH-PFK model variant. The 
model used both generators approximately 28% of the time. They also manifested about 
the same range of use observed in experiment I. We ean eonelude that the model is doing 
what is expected by using smaller generators first to take advantage of better load factor 
for fuel efficiency, which is similar to the results of the first experiment. Generators 1,3, 
and 4, with power ratings of 10 kW, 60 kW, and 10 kW, respectively, are the most utilized 
by the RH-U model variant. All three generators ran approximately 30% of the time. On 
the other hand, generator 2 with a power rating of 30 kW ran the least for both model 
variants, or approximately 16% of the time. This result is consistent with our findings from 
experiment I, and we conclude that this generator’s use requires more investigation as to 
whether a smaller generator with a better load factor could obtain better fuel efficiency. 
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Figure 4.12. Experiment II Factors’ Influence on Cumulative Total Fuel 
Partition Tree, RH-U. 


Proportion of Time Generator Operates 


RH PFK 

RH-U 

Generator 

Proportion in ®/o 

Range 

Proportion in % 

Range 

1 (lOkWO 

28.01 

12.29 -43.19 

29.99 

15.95 -44.52 

2 (30k\\0 

16.26 

3.99 - 36.54 

15.95 

5.32 - 30.23 

3 (60k\\0 

22.50 

0.00-42.19 

30.75 

1.99 - 55.15 

4 (lOkWO 

27.86 

15.61-45.18 

30.03 

17.28 - 45.85 


Table 4.4. Experiment II Proportion of Time Each Generator Operates, 


RH-PFK vs. RH-U. 


Similarly, Table 4.5 reveals that the model rarely requires more than two generators to be 
turned on at any point. Approximately one generator operates 62% of the time for RH-PFK 
and approximately 66% of the time for RH-U. The results show that the power grid almost 
never requires simultaneous use of three or four generators at any point, whieh validates our 
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Effect Summary 
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<.0001* 


Figure 4.13. Experiment II Factors’ Influence on Cumulative Total Fuel 
Stepwise Regression, RH-PFK. 


earlier finding in experiment I. The power grid may operate efheiently with fewer generators. 
This suggestion requires further studies to eonfirm the feasibility in support of fuel savings. 
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Effect Summary 
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Figure 4.14. Experiment II Factors’ Influence on Cumulative Total Fuel 
Stepwise Linear Regression, RH-U. 


Lastly, Figures 4.15 and 4.16 demonstrate our analysis of the load faetor performanee of the 
generators. The quantile summary statistie reveals that 75% of generators 1 and 4’s loads 
have at least an 81% load faetor for the RH-PFK whieh indieates that, when used, these 


51 


















Proportion Number of Generators Turned On 


RH-PFK 

RH-U 

No of Generator 

Proportion in % 

Range 

Proportion in % 

Range 

0 

21.79 

5.65 - 47.84 

13.82 

1.99 - 33.22 

1 

62.10 

43.85 - 76.74 

66.26 

52.82 - 75.42 

2 

15.78 

2.99 - 33.89 

19.08 

8.97 - 32.89 

3 

0.32 

0.00 - 2.66 

0.60 

0.00 - 3.32 

4 

0.004 

0.00 - 0.33 

0.008 

0.00 - 0.66 


Table 4.5. Experiment II Proportion of The Number of Generators Turned 
On At Every Time-Step, RH-PFK vs. RH-U. 


two generators operate at a high load faetor that reduees fuel eonsumption. The situation is 
similar to the RH-U when 50% of generators 1 and 4’s loads have at least a 75% load faetor. 
Generators 1 and 4 run the most for RH-PFK, and they have the best load faetor as well, 
whieh is good for the power grid fuel eonsumption savings. However, the situation differs 
for the RH-U when generators 1, 3, and 4 display almost the same amount of usage, but 
with only 25% of generator 3 load faetor above a 72% threshold. This differenee indieates 
that the generator 3 load faetor is not as good as generators 1 and 4’s load faetors, despite 
being among the most utilized on the power grid. Further study must investigate the power 
grid generator ratings mixture that is best for fuel eonsumption. Our eurrent experimental 
design is not set up to vary generator rating eombinations. 

4.2.5 Other Model Insights 

We now analyze the model temperature behavior within the shelters for additional insights. 
Individual thermostats eontrol the shelters’ temperatures, and the optimization model keeps 
all the shelters within the specified limit except for the four cases in which RH-PFK fails 
to find a feasible solution. Ordinarily, we expect the temperature variation to be narrow or 
cluster since the model’s objective is to minimize fuel usage by keeping the temperature 
within limits. However, the thermostat in shelter 9 keeps the temperature much lower 
than those of other shelters, and shelters 5 and 6 are kept much warmer than the others. 
Figures 4.17 and 4.18 display smoothed results. Smoothing the temperature across 164 
excursions for RH-PFK, and 161 excursions for RH-U, makes overall trends visible that are 
not necessarily evident from the raw data. For instance, we clearly see that shelters 1-4 and 
7 are similar and slightly warmer than shelter 8. A closer look at the raw data reveals that 
shelter 9 has the lowest temperature in approximately half of the time steps across the entire 
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Figure 4.15. Experiment II Distribution of Generators’ Load during Opera¬ 
tion, RH-PFK. 
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Figure 4.16. Experiment 
tion, RH-U. 


Distribution of Generators’ Load during Opera- 
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experiment, and has the lowest average temperature in every single exeursion. We observe 
this behavior in both model variants as well as in both two experiments. As Sprague notes, 
the optimization model frequently maintains shelters below their maximum temperature 
thresholds for long periods of time. This oecurs when a generator must run in order to eool 
other shelters, but its entire eapaeity is not needed for those shelters. Rather than run the 
generator at a low (ineffieient) load faetor for a long period of time, the model ehooses to 
run the generator at a high load faetor for a short period of time, essentially “getting ahead” 
in terms of eooling one of the shelters. Sueh behavior is a natural outeome of a preseriptive 
optimization approaeh, but ean be difficult to encode in a simple policy such as the LGCY 
model. 

One aspect Sprague does not mention is the apparent relationship between the shelter 
characteristics and the average temperature. He states that "the optimally scheduled ECUs 
elects to operate exclusively near the upper, more fuel efficient, limit." However, our results 
indicate that this phenomenon may be at least somewhat dependent on shelter characteristics, 
since different shelters exhibit different temperature behaviors. For instance, shelters 8 and 
9 have different thermal intercept characteristics than the rest of the shelters and tend to 
be cooler, while shelters 5 and 6 have lower-powered ECUs and tend to be hotter than the 
other shelters. Further research could investigate the cause and degree of consistency of 
this behavior among varying grid setups. 

Also, we compare the behavior of the cumulative total fuel consumption for the two model 
variants and the two experiments. This plot represents the averages of the 164 completed 
runs from experiment I, and the 161 completed runs from experiment II. The plot of 
cumulative fuel consumption over time shows a linear relationship with a steeper slope for 
both model variants in the first experiment. However, the relationship is nonlinear in the 
second experiment, with the slope decreasing between time 100 and time 200. This change 
corresponds to the period of day with higher renewable energy production — peaking 
between time step 120 and time step 180, depending on the design point — therefore, lower 
generator usage and subsequently lower fuel consumption. Figure 4.19 displays the results. 


In addition, we examine the fuel consumption variation for five replications across the 33 
design points in experiment II. Figure 4.20 reflects fuel consumption variation for different 
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Figure 4.17. Experiment I Smoothed Temperature Variation, RH-PFK and 
RH-U. 


times of day and different days of the week as represented by the faetor eombinations of 
the design points. This figure helps identify those times or days, as well as faetor variation 
eombinations, responsible for the spike in fuel eonsumption so that decision makers can 
better manage resource allocation accordingly. This insight represents another advantage 
of employing a structured DOE approach with multiple replications across multiple design 
points to gain insight. 
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Figure 4.18. Experiment II Smoothed Temperature Variation, RH-PFK and 
RH-U. 


Figure 4.21 displays the mean relative effieieney of five replieations and 33 design points 
on the left, and on the right shows the same display broken down by design point. The 
data for this eomparison eomes from experiment II, as it uses a more realistie refieetion 
of aetual demand and produetion patterns. This figure presents the times and days as well 
as faetor eombinations of when one model variant performs better than the other on fuel 
eonsumption. Although RH-PFK performs better on average aeross all design points, in 
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Figure 4.19. Cumulative Total Fuel Consumption, RH-PFK vs. RH-U 


some factor combinations on certain times or days, the RH-U performs better. This can be 
due to the particular values of the random data (i.e., “luck”), or an artifact of the rolling 
horizon implementation, or due to the optimality tolerances used in the model; in any case, 
the absolute differences in performance are small in early time steps. However, as the 
time horizon progresses, RH-PFK starts to do better, as expected. A closer look at the 
plot on the right shows design point 6 exhibits wider variation than the others. This is 
the same design point that failed to solve to completion in some instances, indicating that 
it represents a particularly challenging problem instance. Therefore, uncertainty of future 
demand, renewable production, or heat transfer increases fuel consumption, but not by a 
significant margin, as the relative efficiency plot shows. How much a power buffer must 
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Figure 4.20. Cumulative Total Fuel Consumption Variation by Design Point, 
RH-U 


be present for the grid to withstand these uncertainties and improve the power grid system 
resilience is a trade-off decision to be made. 



Figure 4.21. Relative Efficiency, RH-U 


4.3 Analysis Summary 

We have investigated uncertain factors influencing cumulative total fuel consumption by 
analyzing the DOE of a power grid optimization model. Introducing realistic energy 
production and unscheduled demand data into experiment II uncovers additional insights 
that apply in real life situations. Also, our analysis provides useful insights on the generators’ 
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behavior, the temperature’s behavior, and the cumulative total fuel consumption in the 
model. 
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CHAPTER 5: 

Conclusions and Future Work 


5.1 Conclusion 

This thesis investigates uneertain faetors influeneing optimal seheduling of time-shiftable 
electrie loads under uncertainty to gain insights and to support the future development of a 
robust optimization for scheduling such loads. Our analysis of Sprague’s model [1] reveals 
the uncertain factors contributing to fuel consumption, as well as provides useful insights 
on the generators’ behavior in the model. We identify the heat transfer intercept mean of the 
expeditionary shelters as the dominant factor contributing to fuel consumption in the model, 
among the parameters we consider and the magnitude of the uncertainty we investigate. The 
implication is that structures’ or shelters’ characteristics play an important role in how much 
cooling must maintain the shelter’s internal temperature within acceptable limits. Greater 
cooling requirements directly translate into greater fuel consumption. Therefore, choosing 
the right types of shelter when establishing FOBs is a good way to reduce fuel consumption. 

The composition of the energy grid power production system affects the results of the 
identified factor(s) influencing generator fuel consumption. The majority of the power 
production in Sprague’s model is from the generator, as employed in current practice within 
the DOD. However, the outcome varies when we increase the renewable power production 
and unscheduled demand, and then vary them with other uncertain parameters in our DOE 
analysis. Our implementation of these changes in experiment II produces different factors 
that influence fuel consumption. This result implies decisions on setting up FOBs must be 
preceded by a careful analysis and determination of optimal energy production composition 
for the power grid, as well as a realistic estimation of the unscheduled demand requirements 
and renewable energy production. Such practices ensure the right mix of energy production 
capacity prior to establishing FOBs so as to minimize power grid fuel consumption. 

Resource management of limited DOD assets could be improved through insights from 
the analysis of the generators’ behavior. For example, some generators ran at low load 
factors much of the time, even when they were optimally scheduled to run. Although power 
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demand requirements eould increase with little to no notice, our analysis shows improvement 
possibilities, especially with the knowledge of future requirements. Both model variants 
and both experiments illustrate that the power grid requires one or two generators most 
of the time and rarely requires all four generators at once. Running more generators than 
necessary often leads to a low load factor that reduces fuel efficiency. We conclude that a 
structured DOE that investigates power grid generators’ rating combinations, in addition to 
the factors in this study, is required to ensure the right mix of generators is chosen for the 
power grid. 

Analysis of the temperature variations within the expeditionary shelters and the cumulative 
total fuel consumption yields additional insights. The cumulative total fuel consumption 
analysis suggests that uncertainty in unscheduled demand might not lead to a huge drop 
in relative efficiency afterwards. Decision makers need to examine how much to invest 
in a power production buffer for the power grid to deal with uncertainty and increase the 
system’s resilience. Our analysis in this study reveals uncertainty does not harm the power 
grid too much. 

While optimal scheduling of time-shiftable electrical loads has already been shown to reduce 
generator fuel consumption, our analysis of Sprague’s optimization model provides more 
insights on additional improvements. These opportunities could be pursued in different 
ways, ranging from improving structures and shelter characteristics, to changing power 
production composition on the grid. 

5.2 Future Work Recommendations 

Our study reveals fuel consumption of conventional fossil fuel generators could be reduced 
with proper planning and analysis of optimal scheduling of time-shiftable loads. Opportu¬ 
nities exist to improve and continue this research work. 

5.2.1 Robust Optimization Model Development 

Implementation of robustness into Sprague’s model could further reduce generators’ fuel 
consumption, and should be investigated. We performed a preliminary study of developing 
a robust version of Sprague’s optimization model and found that the computation time for 
a robust model is quite high. Future work could investigate this further. 
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As mentioned earlier, some of our findings may support the development of a robust 
optimization model by focusing attention on the aspects of the model and environment that 
make a difference to the solution quality under uncertainty, and the solution timeliness. 

5.2.2 Robustness Investigation via Designed Experiments 

We obtained useful insights from our DOE and an analysis of Sprague’s model. The 
inclusion of experiment II in our analysis offered more insights than we could have obtained 
with only the baseline input data. Other factors, including the power generation capacity 
mix, temperature limit ranges, and so on could be varied to obtain more insights. Future 
work could select other factors to vary, as well as study the behavior of other responses that 
may be of interest for practical applications. 

Also, our analysis provides insights into the behavior of the model; we recommend future 
work directly compare a robust optimization formulation of the model to the DOE approach 
for assessing robustness. 

5.3 Final Thoughts 

This study comes from a desire to improve U.S. military energy use on the battlefield. The 
cost associated with too much energy use goes beyond monetary resources; lives are lost, 
and the budget for future capabilities is diminished. This thesis provides useful insights 
on how to reduce power grid energy consumption by using a powerful combination of an 
optimization model and DOE analysis. Ultimately, the U.S. military expeditionary power 
grid system should be resilient and self-sustaining for our assigned missions around the 
world. 
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